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We provide complete analytic expressions for the inclusive prompt photon 
production cross section in hadronic final states of e + e~ annihilation reac- 



tions through one-loop order in quantum chromo dynamics perturbation the- 



ory. Computed explicitly are direct photon production through first order 
in the electromagnetic strength a em and the quark-to-photon and gluon-to- 
photon fragmentation contributions through first order in the strong coupling 
a s . The full angular dependence of the cross sections is displayed, separated 
into transverse (1 + cos 2 7 ) and longitudinal (sin 2 # 7 ) components, where # 7 
specifies the direction of the photon with respect to the e + e~ collision axis. 
We discuss extraction of fragmentation functions from e + e~ data. 
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I. INTRODUCTION 



Production of an energetic photon in association with hadrons probes the short- distance 
dynamics of electron-positron, hadron-hadron, and lepton-hadron reactions. In addition 
to providing valuable tests of perturbative quantum chromodynamics (pQCD), data from 
electron-positron annihilation reactions permit measurements of parton-to-photon fragmen- 
tation functions. 

In QCD, the quark-photon collinear singularities that arise in each order of perturbation 
theory, associated with the hadronic component of the photon, are subtracted and absorbed 
into quark-to-photon and gluon-to-photon fragmentation functions, in accord with the fac- 
torization theorem Q. Fragmentation functions, D(z,fi 2 ), are inherently nonperturbative 
quantities whose magnitude and dependence on fractional momentum z must be measured 
in experiments at a reference fragmentation scale /ig. The change of D (z,fj, 2 ) with fj, 2 for 
large /i 2 is specified by perturbative QCD evolution equations 0. In e + e~ — > jX, the 
fragmentation contributions play a significantly greater role than they do in hadron-hadron 
collisions || . In lowest-order, the quark-to-photon and anti-quark-to-photon fragmentation 
processes dominate the inclusive reaction e + e~ — > jX, whereas "direct" processes, such as 
qg — ► 7g and qq — > 7*7, dominate in pp — > 7X and pp — > jX for 7's that carry large values of 
transverse momentum 0-0. The dominant role of fragmentation contributions makes the 
inclusive process e + e~ — > 7X a potentially ideal source of information on D (2,/i 2 ). 

The most straightforward theoretical calculations in perturbative QCD are those for the 
inclusive yield of energetic photons, E^da/cPp^. However, an important practical limitation 
of high energy investigations is that photons are observed and their cross sections are mea- 
sured reliably only when the photons are relatively isolated, separated to some extent in 
phase space from accompanying hadrons. Since fragmentation is a process in which photons 
are part of quark, anti-quark, or gluon "jets", it is evident that photon isolation reduces the 
contribution from fragmentation terms. In this paper, we deal with cross sections for ener- 
getic inclusive photons, and the extraction of the photon fragmentation functions D(z,/j?). 
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In another paper ||, we will present a systematic and analytic treatment of cross sections 
for isolated photons, and show over what regions of z the functions D(z,[i 2 ) may be de- 
termined from data on isolated photon production in e + e~ — > 7X. We will also point out 
the breakdown of factorization of the cross section for isolated photons in a particular part 
of phase space in e + e~ — ► 7X, and the implications of this breakdown for calculations of 
isolated photon production in hadronic collisions. 

Our calculations of the inclusive photon yields in e + e~ — > 7X are carried out through 
one-loop order. We compute explicitly direct photon production through first order in 
the electromagnetic coupling strength, a em , and the quark-to-photon and gluon-to-photon 
fragmentation contributions through first order in the strong coupling strength a s . We 
display the full angular dependence of the cross sections, separated into longitudinal sin 2 # 7 
and transverse components (1 + cos 2 # 7 ), where 9 1 is the direction of the 7 with respect 
to the e + e~ collision axis. Our work goes beyond that of previous authors p|-|lT|. For 



example, the full angular dependence of the cross section was not derived before. In one 
recent analysis || , the authors concentrate on events having the topology of a photon plus 1 
hadronic jet; they discuss the extraction of the quark to photon fragmentation function from 
such data. In that approach, final state partons are treated as resolved, and the cancellation 
of infrared singularities is not treated explicitly. Practical aspects of confronting theoretical 



calculations with data from LEP are addressed in Ref. [pLl| . All four groups at LEP have 
published papers on prompt photon production In this paper, we advocate a different 
method of analysis of data from that used so far. 

We begin in Section II with definitions of the factorized inclusive photon cross sections, 
and, to establish notation, we derive explicit expressions for the inclusive photon yields 
in lowest order (O (a° m ) , O (a°)). In Section III, we examine in turn the three 1st order 
contributions to the inclusive photon yield: the O (a em ) process in which a photon is radiated 
from a final quark or antiquark line, e + e~ — > qqj, and the O (a s ) processes in which e + e~ — ► 
qq g, followed by fragmentation of one of the three final-state partons into a photon. The 
quark to photon collinear singularity in e + e~ — > qq^y is absorbed into the quark-to-photon 



fragmentation function. Our treatment of the O (a s ) contributions necessarily includes a full 
discussion of both real and virtual diagrams. Dimensional regularization is used to handle 
infrared and collinear singularities. In Section IV, we summarize our final expressions for 
the inclusive photon cross section E 1 da e + e -^x/d 3 £, with full 6 1 dependence. Numerical 
results and suggestions for comparisons with e + e~ data at LEP, SLAC/SLC, TRISTAN, 
and CESR/CLEO energies are also collected in Section IV. An Appendix is included in 
which we derive expressions for two- and three- particle phase space in n dimensions. 

II. DEFINITIONS, NOTATION, AND LOWEST ORDER CONTRIBUTION 

In this section we establish the notation to be used throughout the paper and present 
our derivation of the lowest order O (a° m a°) contribution to the inclusive energetic photon 
yield in e + e~ — ► jX. 

A. General Structure of the Cross Section and Kinematics 

In e + e~ — > cX, as sketched in Fig. [I], the cross section for an m parton final state is 

dPS (m) -dzD c ^(z), (1) 

with c = "f,q,q,g and z = E^/E c . For inclusive photons, we integrate over all phase 
space, dPS( m \ except the momentum of parton "c". For isolated photons, however, the 
phase space, dPS^ , will have extra constraints due to the definition of the isolated photon 
events. 

For the scattering amplitude, M e + e -^ c+ ..., the vertex between the intermediate vector 
boson and the initial/final fermion pair is expressed as ie^^iyf + a/ 75). The absolute 
square of the matrix element |M| 2 , averaged over initial spins and summed over final spins 
and colors, may be expressed in terms of leptonic and hadronic tensors, and H^, as 

|M| 2 = e 2 C 

4 



da 1 ™' = — 

9e 



PC 



PV 



H 



(2) 



e denotes the electric charge, and C is the overall color factor. Since the physical observable, 
the energetic photon 7, does not distinguish between quarks and antiquarks, the parity 
violating (PV) term does not contribute. Equivalently, only the symmetric part of H^ v 
contributes. Therefore, 

|M| 2 = e 2 CF PC (q 2 ) = e 2 CF q PC (q 2 ) {H l + H 2 ) . (3) 

Hi = {-9,u + W» = -g, u H^. (4) 

H 2 = -^H^. (5) 
q z 

The four-momenta g M and k^ are defined in terms of the four-momenta of the incident e + 
and e~ (fcfand k 2 ) as 

q» = k? + k$, q 2 = (h + k 2 f = s; (6) 

and 

jfe" = A£ - k%, k 2 = (ki - k 2 f = -s. (7) 

The normalization factor F PC (q 2 ) is expressed in terms of the vector (v) and axial- 
vector (a) couplings of the intermediate 7* and Z° to the leptons and quarks. At the Z° 
pole, neglecting 7, Z° interference, we find 

\ Fro - (ki- + w) h + w) (s _ MI ; +MIrr («) 

At modest energies where only the 7* intermediate state is relevant, 



- (s) = e\ i; (9) 
s s 



e q is the fractional quark charge [e u = 2/3; = 1/3; ■ ■ •). 

In terms of functions Hi and H 2 , defined through Eq. (JJ), we reexpress the cross section 



as 



da™ = £ 



-F, 



PC I 



e 2 C - (Ht + H 2 ) dPS {m) dzD(z) 



(10) 



In the sections to follow, we calculate the functions Hi and H% explicitly for the lowest 
order and the first-order contributions to e + e~ — > jX. Function Hi provides the cross 
section integrated over all production angles, 9 y , of the 7. Function H 2 specifies the angular 
dependence or, equivalently, the transverse momentum distribution of the 7 with respect to 
the e + e~ collision axis. 



B. Factorized Cross Section 

We are interested in the inclusive cross section for production of photons in association 
with hadrons, E y da l Jl c ^^ x /d 3 £, where E y is the energy of the photon, and i is the momen- 
tum of the photon in the e + e~ center-of-mass system. According to the pQCD factorization 
theorem M, we may express the cross section as 

e£^^E*^*D^). (11) 

The intermediate partons are c = j,g,q, and q. The hard-scattering cross section 
E c da™ c ^^ cX /d 3 p c contains no infrared or collinear divergences. The fractional momen- 
tum z is defined as z = E 1 /E c ] all intermediate partons c are assumed to be massless. 
The fragmentation functions D c ^y(z) represent all long-distance physics associated with 
the hadronic component of the photon. They are inherently non-perturbative quantities 
that must be measured experimentally. Models and phenomenological parametrizations WA 



for D(z) have been published. In lowest-order, Z) 7 ^ 7 (z) = 5(1 — z). The convolution ex- 
pressed in Eq. flTTp is sketched in Fig. |2|. The symbol ® in Eq. fllTD is defined explicitly as 
follows: 



d s p c 



2 



'da^ cX (E c 

E r 



D c ^{z). (12) 



d 3 p c 

Since z min occurs when p c has its maximum value, p™ ax = y/s/2, the lower limit of integration 
Zmm = x,y = 2Ej/ y^; yfs is the center of mass energy of the e + e~ annihilation. 



C. Derivation of the Lowest Order Contribution 



In this section we present an explicit derivation of the lowest order contribution to the 
inclusive photon yield in e + e~ — > 7, sketched in Fig. [5| The differential inclusive cross section 



da P 



+7 x is expressed as a product of the lowest order partonic cross section da^ e ^^ q - 



and the q — » 7 fragmentation function, D 



da, 



e+e~— y^X 



5^7 V 



e+ e — >q(p q )<i dz D P^j( Z ) + (?->?)■ 



(13) 



In Eq. (|T3"D, p q is the four- vector momentum of the quark q, and z = £/p q . The partonic 
cross section is written, in turn, in terms of the invariant matrix element and differential 
phase space factor. 



da 



(o) = J_ 

2 



F„ PC ( 



*PqPq 



dPS® 



s) 



e 2 Ar c -(iy 1 + J H- 2 ) dPS^ 



(14) 



where N c = 3 is the number of colors carried by the quarks, and Eq. (|3]) was used. 

The symmetric part of the hadronic tensor H^, used to define functions Hi and H2, is 
particularly simple: 



= 4 (e//) 2 + P%P« ~ 9" v P q ■ V~ q 



(15) 



The factor // in Eq. ( |I5j ) accommodates the fact that we are working in n dimensions. The 
dimensional scale fi will be specified further below. The functions Hi and if 2, defined in 



Section [II A| , become 



Hi = A{e^f s(l-e); 
H 2 = -2(e//) 2 s(l-cos 2 i 



(16a) 
(16b) 



In Eq. (|15D , e is defined through the number of space dimensions n = 4 — 2e, with e — > at 
the end of the calculation. In the center of mass frame of the collision, 9 is the angle of p q 
with respect to the direction defined by the incident e + . Combining Hi and if 2, we obtain 



~ {Hi + H 2 ) = \ (e//) 2 s 



[1 + cos 2 



2e 



(17) 



Combining Eqs. ( |i~7|) and the expression for two-particle phase space in n-dimensions, 
Eq. ( |A4l ) of the Appendix, we find that the lowest order partonic cross section, Eq. (|T4|), is 



do - + - x- 

9 <P Pq 



2 / 47^ . x 

em c ^/4)sin 2 ^ r(l-e) 



(l + cos 2 9) - 2e 



8(x q - 1) 



X (l 



(18) 



with x q = 2E q / \fs. At this order, the cross section is manifestly finite in the limit e — > 0, 
and we may set e = directly in Eq. (0). Nevertheless, Eq. flT8| ) expressed in n dimensions 
is valuable for later comparison with the higher order cross section. 

Noting that t = zp q implies d 3 p q / E q = (l/z 2 ) d 3 £/ E 7 , we obtain the lowest order inclusive 
cross section 



7 dH 



da 2 - v 

""^ e + e -_ y qX 

q d?p q 



X r , 



2E 



-F q PC (s] 
s H 



z 

1 



a em N c (l + cos 9 1 )—D q ^ 1 {x 11 fi F ). 



(19) 



The angles Q 1 and 9 are identical since we take all products of the fragmentation to be 
collinear. The overall factor of 2 in Eq. (|T9|) accounts for the q contribution. In Eq. ([19]) we 
have introduced a fragmentation scale \ip in the specification of the fragmentation function. 

III. FIRST ORDER CONTRIBUTIONS 

There are three distinct contributions to e + e~ — > jX in first order perturbation theory: 



e + e -> 7, 0(« er 

e + e" -> g (or g) -> 7, 0(a s ) 

and 

e+e" -> g -> 7. 0(a s ) 



(20a) 
(20b) 

(20c) 



S 



In Eqs. (|20b|) and ( |20cj) , we have in mind contributions from quark and gluon fragmentation 
to photons in the three-parton final state process e + e _ — > qqg. The first contribution, 



Eq. ( |20a| ), arises from e + e~ — > qcpj where the 7 is not collinear with either q or q. 

In this section we derive and present the explicit contributions to the inclusive yield 
E^da™^- _^ yX / 'd 3 £ from each of the three processes in Eq. (20). Following the pQCD factor- 
ization theorem, and Eq. (|ll|), we must calculate the short- distance hard-scattering cross 
sections, E c da™ d e -^ cX / 'd 3 p c for c = 7, g, q and q. 

The Feynman graphs for e + e~ — > ^qq are sketched in Fig. Owing to the quark-photon 
collinear divergence, the cross section associated with these graphs is formally divergent. We 
denote this divergent first order cross-section (T^ e - x , a short-hand notation for Eda/d 3 £. 
To derive the corresponding short- distance hard-scattering cross section, <7g+ e - x , we apply 
the factorized form, Eq. ([□]), perturbatively, 

^l-^x = tfl-^x ® (z) 

+ (<Z-g). (21) 

The convolution represented by ® is defined in Eq. flT2l) . The superscripts (0) and (1) on 
the hard-scattering cross sections a and fragmentation functions D refer to lowest-order and 
first order, respectively. The collinear divergence resides in the first order fragmentation 
function D^X^z). The hard-scattering cross sections and are finite. The expression 

7 x : 



for o~(°) was derived in Section |II C| . In Section A we present our derivation of af} i 



^Se-^x = ^e--nx - ^l-^x ® D(%(z) - (q - q). (22) 

The two Feynman graphs that provide the cross section for e + e~ — > g — > 7 in 
0(a s ) are shown in Fig. |5]. In this case, the final gluon is effectively "observed" through 
the fragmentation g —>■ 7; there are no virtual gluon exchange diagrams. The finite hard- 
scattering cross section cr^+ e -^ gX is derived from the difference 



= - t - { l-^x ® D% g . (23) 



In Eq. (p3|) , the divergent cross section cr^+ e -^ gX is evaluated from the Feynman graphs 
shown in Fig. [5], and the quark-to-gluon collinear divergences are embedded in the first- 
order fragmentation function D q ,_^ g . Our derivation of ffg+ e - x is presented in Section B. 

The Feynman graphs in 0(a s ) that contribute to e + e~ <1 ^ 1 (Eq. (|20 b| ) ) 
are sketched in Fig. |6]. Only a final state photon from quark fragmentation is observed. 
The complete 0(a s ) result includes both real gluon emission and virtual gluon exchange 
graphs, as shown in Fig. ^|. Although infrared divergences associated with soft gluons cancel 
between the real and virtual graphs, the cross section <7^+ e - x obtained from the Feynman 
graphs is still divergent due to collinear singularities when the real gluon is emitted along the 
direction of its parent quark or antiquark. To obtain the corresponding hard-scattering cross 
section 6"e+ e --*gX> we apply the factorized form, Eq. (0), perturbatively, to the production 
of a quark instead of the photon, 

° { X-^ X = i eu, ® D% q + eu 9 * ® (24) 



1 —q 



with the collinear q' — > q singularities in 0(a s ) included in DQ' . Note that Dj® 



q'^q- " " q^q 

(1) 



5(1 — z). Correspondingly, the finite hard-scattering cross section (t e + e -^ qX is 

<r (1) -rx (1) -<r (0) (251 

& e + e -^ q X - a e+e~->qX ^e+e"-^ ® U q^q' l Z0 > 

In Section C we present a detailed derivation of cr^+ e - x . 

Before turning to our explicit derivations, we conclude this discussion with a presentation 
of the factorized formula for the two-loop short- distance hard-scattering cross sections. The 
two-loop direct contribution to e + e~ — > 7X is of 0(a em a s ) and can be derived as follows. 
We first apply the factorized form, Eq. (|iTl), perturbatively, at two- loop level and sum over 
c = J,g,q and q, 

u e+e~-^jX 













-^ 7 x c 










» 














» 









+ (26) 
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All first-order contributions, cr^\ in Eq. ( f26|) are given in Eqs. (^), (|23|) and (|25|) , and 
they are calculated in this paper. Since the first order fragmentation functions I)W (z) and 
D^Xj(z) vanish, and the zeroth order hard-scattering cross section &f+ e -_> g x vanishes, we 

(2) 

derive the two-loop hard-scattering cross sections v e + e - _+ lX as: 

~ (2) _ (2) 

u e+e-->jX ~ u e+e-^jX 

-(_-_). (27) 

To complete the calculation of cr e + e -_^ lX , it is necessary to calculate the two-loop parton-level 
cross section cr+ e __> x and the two-loop quark-to-photon fragmentation function D^ 7 (z) 
in n-dimensions (implicitly, we use dimensional regularization), in addition to all the zeroth 
and first order contributions calculated in this paper. The two-loop parton level cross section 



er+ e - x is formally divergent. As is true of the calculation of <7g+ e - x in Section |lll Q , all 



infrared divergences associated with soft gluons cancel among the real emission and virtual 
exchange diagrams. All collinear divergences that appear when final-state quarks and/or 
gluons are parallel to the observed photon are cancelled by the subtraction terms given in 

_____ (o\ 

Eq. fl27|). Consequently, the two-loop hard-scattering cross section a e + e - x is finite if the 
pQCD factorization theorem holds. 

As shown in Section IV, the leading order short-distance direct production con- 
tribution <7g+ e - x is much smaller than the leading order fragmentation contribution 
cfe+ e - x ® D^X^z) + (q <-» q). We expect that the next-to-leading order direct contribution 

A (2) 

Og+g-^jf will be much smaller than the next-to- leading order fragmentation contributions 
^e+e--,cx ®D2Xj( z ) with c = g,q and g, which are completely derived in this paper. We will 
not calculate the two-loop contributions in this paper because we believe their contributions 
to the overall cross section are much too small in comparison with those presented here. 
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A. Derivation of a 



(i) 

e+e — — +7X 



In this section we present an explicit derivation of the finite hard-scattering cross section 
0"e+ e --vyx * n ^W' We begin by computing the functions Hi and H 2 , defined in Sec- 
tion [11 A| , Eqs. M) and (|5|). These will then be integrated over phase space to yield the cross 



section E^d<j^+ e - x / 'd?£ 
da { }} 



E 



-F q P %) 
s y 



e 2 N c h y Hi + H 2 )dPS^\ 



(28) 



where three-particle phase space in n-dimensions is given in Eq. (|A29|) of the Appendix. 

Sketched in Fig. |7| is the hadronic tensor H^ obtained from the two diagrams of Fig. |j. 
Performing traces to sum over final spins, we may write the four contributions as 

1 



H$ = 2(1 -e)Tr [wlwp*] 



2pi 



Hj£ = 2(1 -e)Tr [ 7/l7 • p llvl ■ 



2p 2 



H l C u = ~ 2Tr [7^7 • Pi7 • P2lul ■ (Pi + ih ■ (P2 + 



+ 2eTr [ 7/i7 ■ p l7 ■ £ 7 ^ 7 • p 2 7 ■ 



1 



1 



2pi -£2p 2 - 



2pi ■ £ 2p 



'2 ■ 



HjfJ = -2Tr [ 7 „ 7 • {pi + t)<y ■ (p 2 + £) lvl ■ p l7 • p 2 ] 



+ 2eTr [ 7m7 • £ 7 • p i7l , 7 • £ 7 • p 2 ] 



1 



1 



2pi • £ 2p 2 



(29) 



2px • £ 2p 2 ■ £ 

To obviate multiple repetition of a common factor, we temporarily omit the overall coupling 

d 

factor e;j(e//) 4 that appears in H^ u . Function Hi = —g^ u H^ v = —g^ J2 H { - l ' tMU . We obtain 



Hi = 8(1 - e) (1 - 6) 



Vie + me 

V2t Vu 



+ 



2y 



12 



2e 



VU V2i 

The dimensionless quantities yu, y 2 £, and yi 2 are defined by 

2ft ■£ , 



(30) 



2/12 



2pi ■ p 2 
q 2 



1,2); 



(31) 
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We remark that yu + yu + Vii = 1- In evaluating H 2 = — (k^k u /q 2 ) H^, we also make use 
of dimensionless quantities yik,V2k, and y^f 

2 Pl ■ k . 



2k ■ 

q 2 



1,2 ; 



Because k ■ q = 0, 



Vik + Vik + Vu = 0. 



(32) 



(33) 



After some algebra we find 

H 2 = -A(l-e) 



yu + y2e_ 

V2i Vu 



+ 



2y 



12 



2e 



Vu Vii 



+ 



{ylk + vlk} — — {yu} ■ 

<- > 111 nil nn y J 



Our next task is to integrate H\ and over three-body phase space in n = 4 — 2e 
dimensions. Since the momentum of the photon (£) is an observable, and the momentum of 
either the quark (pi) or antiquark (p 2 ) can be fixed by the overall momentum conservation 
^-function in the three-body phase space, we need to integrate over only pi or p 2 . In the 
following discussion, we let p 2 be fixed by the 5-function, and we integrate over p\. In the 
overall center of mass frame, as sketched in Fig. || we take angle # 7 to be the polar angle of 
the 7 with respect to the e + e~ collision axis and angle #i 7 to be the angle between the 7's 
momentum I and the quark momentum p\. The angle 9 X is the n-dimensional generalization 
of the three-dimensional azimuthal angle 0, defined through p\ as 



dVL n _ 2 (pi) = d9 Xl sin n 3 9^d9 x sin n 4 9 x dVt n _ 4 (pi). 



(35) 



Having chosen the frame, we may reexpress the y variables in terms of observables and 
integration angles as follows: 



2/w2/i2 - yu 



yik 
y-ik 



,1 ~, 



yuyi2 - y2i 



cos 9~ 



cos 9~ 



2 v /2/l22/l£2/2£ 



sin 9 j cos 9 X ; 



sin 9 7 cos 9 X ; 



(36) 
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where x 1 = 2E y /y/s(= yu + V2i)- In deriving y lk and y 2 k, we use the following identities 

y2tyu - yu 



cos#i 7 
sin 6*i 7 



2y/yi2yuU2t 



(37) 



where Xi = 1E\j \fs. 

In the integration of H 2 over phase space, the integral over dcos6 x is done from cos^ = 
— 1 to +1. The expression for the three-body phase space, Eq. ( |A23| ), is an even function 
of cos9 x . Correspondingly, terms in H 2 that are odd functions of cos9 x do not survive. 
Because H 2 depends only on the square of the yik and y 2 k, after eliminating all terms linear 
in cos we find that the only X dependence in H 2 is cos 2 ^. We can integrate over 9 X 
independent of other variables, or we can effectively replace the cos 2 9 X terms in H 2 by the 
average of cos 2 9 X in n-dimensions and eliminate the 9 X dependence in H 2 completely. 

Given the average of cos 2 ^ in n-dimensions, Eq. ( |A26|) , we obtain, effectively, 



2 2 2 n 



7 ' 



ylk 



ylk 



yuy\2 - y2i 



1 2 



COS 2 9y + 



cos 9y + 













^{y\2VuV2t) 




X 2 



sin 2 6L ; 



sin 2 9y ; 



(38) 



where the factor 1/(1 — e) is from the average of cos 2 9 X . Substituting the above expressions 
into Eq. (0), and combining with Hi in Eq. ([$(]), we obtain, 



j (Hi + H e 2 n ) = (l + cos 2 9y - 2e 



{1 _ e) yu + yju 

\y2t yu 



M 



_yi2_ 

yuyn 



+ 1-3 cos 2 9r 



+ 



4 2/i2 



X 



7 J 



cos 9- 



4 2/i2 



X 



7 J 



(39) 



where the superscript "erf' indicates that we have replaced cos 2 9 X by its average in n- 
dimensions. The two ^-functions in the three-particle phase space, dPS^ 3 \ provide the 
following identities 
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yi2 = 1 - x-y ; 
V2t = x 1 - yu ■ 

Introducing yu = yu/x-y, and substituting these identities into Eq. (1351), we derive 



(40) 



#1 + H e 2 ff 



1 + cos - 6L — 2e 



+ (l + cos 2 # 7 - 2e 



+ (l - 3 cos 2 ^ 



1 + (1 -x 7 f 



xi 



-2 - e 



1 



1 



1 



1 ' 



4(1 -x 7 ) 



xi 



(l - cos 2 ft. 



4(1 



(41) 



1-e, 

The last term vanishes as e — > 0. 

Combining Eqs. ( ^8|) and fl4l|), and integrating over <i?/i^, we can derive the partonic cross 
section der+ e _ x . The limits of the dyu integration are from to 1. The integrals over 
dyu for (Hi + H% )/4 may be expressed in terms of 



he Vu * (! - m) 



Examining Eq. (f4l"D, we need only J 0j o and I-i } o (= -^0,-1) : 



(42) 



/ ,o = 5(l-e,l-e) 



J_x n = B(-e, 1 - e) 



_ (r(i-6)) 2 , 

r(2-2e) ' 

1 \ (r(i-e)r 
- e ; rn-2e) 



(43) 



(44) 



For small e, / ,o — 1 + 0(e), and I_i j0 = — \ + 0(e). 

After performing the integration over dyu, we expand the right-hand side of Eq. fl28|) 
in a power series in e, keeping only the singular term proportional to (1/e) and the terms 
independent of e. (Terms of 0(e m ), m > 1, vanish in the physical limit of four dimensions 
(n = 4 - 2e)). We obtain 
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da 



(i) 



d 3 £ 



2E 



2 / 47T/1 2 \ 6 1 1 / 2 



x- 



+ 2E 



2 "em 

q "27 



F„ PC ( 



1 + (1 - 



x|(l + cos 2 7 
+ (l - 3 cos 2 ft 



.1 



1 + (1 - X. 



1 

e 

~2tT 



2(1 -x. 



(45) 



In deriving Eq. (^5p, we included the overall factor for coupling constants, e 2 (e//) 4 ; and 
used the expansion T(l — e) ~ 1 + e^E, where 7^ is Euler's constant, and the usual modified 
minimal subtraction scale 



-7-B 



(46) 



The ( 1/e) singularity in Eq. ( P>D represents the quark-photon collinear singularity. This 
singular term is expected to be cancelled by subtraction terms defined in Eq. fl2"2"|). By evalu- 
ating the diagram sketched in Fig. ^, we obtain the one-loop quark-to-photon fragmentation 
function 



9-*7 



2 "em 

3 " ~27 



i + d 



(47) 



where we keep only the 1/e pole term because we work in the MS factorization scheme. 
Using the fact that D^X^z) = D^l (z), and comparing Eq. (flop with Eqs. (0) and (fl7|), 
we observe that the divergent first term in Eq. (^) is cancelled exactly by the subtraction 
terms defined in Eq. (|22|) , in accord with the pQCD factorization theorem. Using Eq. (p2|). 
we obtain the finite 0(a em ) hard-scattering cross section 



da 



(i) 

e+e — — *-yX 



dH 



2E 



iPC , 



x| (l + cos 2 ft 
+ (1 -3 cos 2 7 ) 



;i-x 7 ) 2 



2(1 - 



'' * 2tt 

£n (s/ /x|jg) + in (x 2 (1 — x n 



(48) 
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We remark that the angular dependence of the 0(a em ) hard-scattering cross section has 
two components, one proportional to (1 + cos 2 9), familiar from the lowest order expression, 
and a second piece proportional to (1 — 3cos 2 9 1 ). If one integrates over cos# 7 , the second 
piece vanishes. We note, however, that the piece proportional to (1 — 3cos 2 6> 7 ) changes 
the predicted angular dependence from the often assumed form (1 + cos 2 9). The difference 
means that it would not be correct to assume a (1 + cos 2 #) dependence when attempting 
to correct an integrated cross section for unobserved events near, for example, the incident 
e + e~ beam direction, cos# 7 = ±1. 

B. Derivation of a + _ v 

e+e — >gX 

The finite hard-scattering cross section 0"e+e--»sX ^° ^ rs ^ or der in a s may be obtained 
directly from Eq. fl4§D after three replacements: x 7 — > x g ; N c — > N c Cp\ and e 2 e 2 of the final 
photon emission vertex by g 2 = 4ira s . 



E, 



2E 



F„ PC (a) 



a 2 em N c 



2tt 



x • 



(l + cos 2 9 t 
1 - 3 cos 2 9, 



l + (l-a; ) 2 



Xg 

2(1 -x. 



X r . 



in (s//i|jg) + in (x 2 g (1 - x g )) 



(49) 



In Eq. (H), x g = 2E g /^/s; C F = §, and N c = 3. 

The contribution 0(a s ) to the inclusive yield e + e~ — > 7X via gluon fragmentation is 
therefore 



E. 



(i) 



e+e -->gX^yX 

dH 



1 dz 



9 <P Pg [ X9 ~ Z 



(50) 



with x 7 = 2E 1 /y/s. Because the g — > 7 fragmentation process is collinear, 9 g = 9 T 



C. Derivation of a 



(i) 

e+e-^qX 



In this section we present our explicit derivation of the finite hard-scattering cross section 
^e+e--*qx/ c ^Pq to first order in a s . As sketched in Fig. || both real gluon emission and 

17 



virtual gluon exchange graphs contribute. The real emission diagrams have both infrared 
and collinear divergences. The infrared divergence is cancelled by contributions from the 
virtual diagrams, while the collinear divergence is cancelled by the subtraction term defined 
in Eq. 

The real emission diagrams can be treated easily in the same way as da^+ e _^^ x /d 3 £ 
in Section |III A| . Except for the replacement of a photon by a gluon, the hadronic tensor 
H^y obtained from the gluon emission diagrams in Fig. ^a is identical to that computed in 



Section [III A| for e + e~ — *> qq^y. Thus, we may employ our previous expressions for Hi and 
H 2 again but with the replacement of subscript "f in Eqs. ( |30D and (|34|) by subscript "3" , 
since p^ is our momentum label for the gluon. Because the quark is now the fragmenting 
particle (i.e., effectively the "observed" particle), the yf k variables with i = 1,2, 3 in H 2 are 
no longer those in Eq. (0). Instead, we now have 



z 2 2/i 



1/131/23 - 1/12 



X 1 



and 



1/121/23 - 1/13 



Xi 



cos 2 9i + 



cos 2 6>i + 



1 -e 



1 - e 



2(1/121/131/23) 



2(1/121/131/23) 



sin 2 a 



sin 2 9 1 



(51) 



In Eq. (^Tj), 6>i is the scattering angle of the quark, and subscript "3" indicates the gluon 
of momentum p^. We dropped all terms linear in cos^, and replaced cos 2 9 X by its average 
value in n-dimensions. Substituting these yf k with i = 1, 2, 3 into Eq. fl34D , we derive 



eff 
2 



1 + cos 2 - 2e 



(1 _ e) yjz + y™ 

Vl/23 1/13 



M 



+ 1 - 3 cos 2 e 1 



21/12 



x\ 



+ e cos 2 #1 



4 1/12 
x 2 



where the last term again vanishes as e — >• 0. In analogy to Eq. 
here are 



I/12 

1/131/23 



(52) 

the useful identities 
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2/12 = Xi - J/13 



(53) 



Using these identities, we reexpress Eq. 
1 



4 



1 + cos^ Q x - 2e 



in terms of X\ and 1/13 

1 + 2? 



+ (l - 3 cos 2 0i 



Xi 



1 — Xi 



yi3 



1 

2/13 



2/13 



\ — X\ 1 — X\ 



2/13 



1 — Xi 



• r l 



(54) 



where we dropped the last term in Eq. (52). Introducing the overall coupling factor 
(e/i e ) 2 (g//) 2 and color factor N c Cp, and combining with the three particle final state phase 
space dPS^ 3 \ Eq. ( |A3C ), we express the contribution of real gluon emission as 



£1 



d 3 pi 



x CV 
1 

x - 
4 



2tt 



a 2 AT 



47T/i 2 



'47r/i 2 



(s/4)sin 2 0iy r(l-e) 

5 (xi - (1 -2/23)) 



s J r(i 

^2/i2 c?2/i3 dy: 



• J 'i 
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^ (1 - 2/12 - 2/13 - 2/23) , 



(55) 



2/12 2/13 2/23 

where superscript (R) stands for the real emission. Using the two (^-functions to fix j/23 and 
2/12, inserting Hi + B.^ from Eq. (|54|), and integrating dy\% from to xi, we derive 



da 



(R) 
e+e" 



d 3 pi 



?PCi 



a 2 N 

"em 1 v c 



47T/i 2 



X 



'47T/i 5 

1 + cos 2 0i - 2e 

/l +x 2 \ 
+ ' 1 



^/4)sin 2 ly / T(l-e) 



T(l-e) 



r(l-e)J Vxi/ r(l-2e) 
7 ^ r + 3 -S(l-x 1 )" 

l-X 1 ) + 2 y 

2 \ / £n(l — xi) 



£n (xf + 1 + xf 



1 — Xi 



1 — Xi 



2 VI - xi 



1 -3 cos 2 0x) 



(56) 
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The "+" prescription is defined as usual 



(57) 



The right-hand-side of Eq. (|56| ) is formally divergent as e — > 0. The 1/e poles in n-dimensions 
represent the infrared divergence, when the gluon momentum goes to zero, and/or a collinear 
divergence, when the gluon momentum is parallel to that of the fragmenting quark. As we 
show below, the infrared divergence is cancelled by the infrared divergence of the virtual 
diagrams, sketched in Fig. ||b. 

The contribution of the virtual diagrams results from the interference of the one-loop 
vertex and self-energy diagrams with the leading order tree diagram. As for the leading 
order contribution, the virtual diagrams, sketched in Fig. ^b, have a two-particle final state 
phase space. Therefore, the contribution from the virtual diagrams has the same kinematical 
structure and angular dependence as the leading order contribution, discussed in Section |I C| . 
It is proportional to 5(1 — x±), and, consequently, the virtual contribution cancels only the 
1/e poles associated with the 5(1 — Xi) terms in Eq. (|56|). The subtraction terms in Eq. ([25]) 
cancel the final state collinear poles that appear in the contribution of real gluon emission. 

Beginning with the virtual exchange diagrams in Fig. ^|b, we evaluate the one-loop vertex 
correction in n-dimension, and combine it with the lowest order tree diagram to form the 
first order virtual contribution. We derive 



E 1 - 



do 



(V) 

e+e-^qX 



d 3 pi 



-F PC (s) 
s 



O's 



a 2 N 



(s/A)sm 2 e 1 ) r(l-e) 



47T/1 



2\ e 



r(i-e) 



1 + cos 2 9 X -2e) 5(1 -xx 



i_\ r(i-e) 3 r(i + 6) 

xj r(l-2e) 
2 3 / 2 

-— h 7T - 



(58) 

where the superscript (V) stands for the virtual contribution. After adding the real and 
virtual contributions, Eqs. (0) and (|58|), we obtain the cross section for e + e~ — > qX at 
order 0(a s ), 
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4-SL 



s 



a "" JVc l( S /4)sin 2 eJ r(l-e) 



.PC, 



a em iV c 



an 



2tt 



x<^ (l + cos 2 0/ 



J±|L + 3, (1 _ Xl) W_£ 



MS' 



2 VI - an 



+5(1 -xi] 



2tt 2 9 



(3xi — 5) 



+ (l - 3 cos 2 0i) | 



(59) 



As is evident from the 1/e terms, this cross section is divergent as e — > 0, a reflection of the 
fact that a cross section for producing a massless quark is an infrared sensitive quantity, not 
perturbatively calculable. 

According to the pQCD factorization theorem, the short-distance hard-scattering cross 
sections, defined in Eq. ([11]), are infrared safe quantities. Beyond the Born level, the short- 
distance parts, a e + e -^ c x, are not the same as the partonic cross sections a e + e -^ c x for frag- 
menting parton c. Following Eq. (p5j), in order to derive the short-distance hard-scattering 
cross section o"e+ e --+gX> we must first calculate the one- loop perturbative fragmentation 
function . Feynman diagrams for D^\ q are sketched in Fig. 10. These diagrams are 
evaluated in the same way as one evaluates parton- level parton distributions [[0|, and we 
obtain 



2tt 



l + x\ 3 r/1 

;T^T7 + 2 5(1 -^ 



(60) 



where the "+" prescription is defined in Eq. (|57|). 

Using Eq. fl25|) , the lowest order cross section for e + e~ — > qX, Eq. flT8|), and the one-loop 
quark fragmentation function D^\ ql Eq. (|BU[), we derive the short-distance hard-scattering 
cross section 
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<»S> 



d 3 pi 



a 2 N — 

171 x 



x (1 + cos 2 6- 



in 



+ 



'1 + gf' 

1 — X\ 

-5(1 -x x ) 



in (x\) + (l + xl) 



in(l — Xi) 



1 — Xi 



2 VI -xi 



+ (l -3cos 2 # 7 ) |. 



(61) 



We set # 7 = $i based on the assumption of collinear fragmentation from quark to photon. 
As expected, the hard-scattering cross section is infrared insensitive. The 0(a s ) quark 
fragmentation contribution to e + e~ — > 7X is 



d 3 £ 



q J X'y Z 



£1 



d 3 pi 



x, 

Xi = — 

z 



(62) 



Our derivation shows that the short-distance hard-scattering cross section for antiquark 
fragmentation to a photon is the same as that for quark fragmentation. Consequently, the 
0(a s ) antiquark fragmentation contribution to e + e _ — > 7X is the same as that given in 
Eq. (M). 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we present and discuss explicit numerical evaluations of the inclusive 
prompt photon cross sections derived in this paper. We provide results at e + e~ center-of- 
mass energies y/s = 10 GeV, 58 GeV, and 91 GeV appropriate for experimental investigations 
underway at Cornell, KEK, SLAC, and CERN. In our figures, we display the variation of 
the inclusive yield with photon energy E~ and scattering angle # 7 , where 9 1 is the angle of 
the photon with respect to the e + e~ collision axis. We also show the dependence of cross 
sections on the choice of renormalization scale \i. 

The cross sections we evaluate are those derived in the text: Eqs. (0), (f%8"D, Q5TD, and 



). They are assembled here for convenience of comparison. The lowest order inclusive 
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cross section is 
daiTi 



E. 



7 



e+e ^yX 



2E 



f: c {s) 



U 2 em(s)N c (l + COS 2 OJ—D^X^ fij,). 



The finite 0(a em ) hard-scattering cross section is 
7 ~ dH 



do + 

77, e+e ~^yX 
/'.- 



2E 



2 ( a em{l^F) 



(63) 



2tt 



x | (l + cos 2 7 ) 
+ (1 -3cos 2 # 7 ) 



1 + (1 - x 



2(1 -Xy) 



in (s/fi 2 F ) + £n (a; 2 (1 - x 7 )) 



(64) 



The 0(a s ) contribution to the inclusive yield e + e — > 7X via gluon fragmentation is 



» e+e-^gX^ 7 X 



1 (fc 



9 <Pp g 



X n 



D g ^ 7 (z,n 2 F ) 



(65) 



with 



E n 
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e+e-^gX 
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x I (l + cos 2 6> 7 ) 
+ (f - 3 cos 2 7 ) 
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2(1 -x fl 



'a s (/x 



2tt 



£n (s//4) + in(x 2 g (l-x g )) 



(66) 



We choose the renormalization scale \i in a s (/i 2 ) to be the same as the fragmentation scale 
Hp in D g ^(z, up). The 0(a s ) contribution to the inclusive yield e + e~ — > 7X via quark 
fragmentation is 



da 



(i) 

e+e-->gX-> 7 X 

dH 



1,1 cte 

q J Xj Z 



E d °e+e-^ q X 

d 3 p! 



D q ^(z,^i 2 F ) 



(67) 
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with 



da 
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e+e-^qX 



d 3 pi 



a lm( S ) N c — 



c f 



2tt 



x (1 + cos 2 # 7 ) 



+ 



+6(1 -xi) 



1 - 3 COS 2 ft. 



2tt 2 9 



27-2 (3xi - 5) 



(68) 



For the common overall normalization function F^°(s), we use an expression that in- 
cludes 7, Z° interference: 



4+ (\ l 'e\ 2 + K\ 2 ) (\Vq\ 2 + 



(s - Ml) 2 + Mfr| 



s ( s -m 2 ; 



" ' "s-M%) 2 + MlY\ 



(69) 



The vector (t>) and axial-vector (a) couplings are provided in Table | and Table [TT]. We set 
Mz = 91.187 GeV and Tz = 2.491 GeV. These and other constants used here are taken 
from Ref. ||15|| . The weak mixing angle sin 2 6 W = 0.2319. For the electromagnetic coupling 



strength a em , we use the solution of the first order QED renormalization group equation 



Qiem(/^ 



1 + ga em (/i 2 )^(/xV/i 2 ) 



(70) 



Here /?o is the first order QED beta function 
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(71) 



/ 



with N[ the number of colors for flavor / and ej the fractional charge of the fermions. The 
sum over / extends over all fermions (leptons and quarks) with mass m 2 < /i 2 . For the 
energy region of interest here, we do not include the top quark in the sum in Eq. ([71]), and 
we obtain /3 = —80/9. To fix the boundary condition in Eq. ([T0|), we let a em (M|) = 1/128 
and set /xq = Mz- 



24 



In the O (a s ) contributions, Eqs. (|66D and (0), we employ a two-loop expression for 
a s (/i 2 ) with quark threshold effects handled properly. We set Ag^ D = 0.231 GeV. At 
y's = Mzi this expression provides a s (M§) = 0.112. 

At yfs = 10 GeV, the sums in Eqs. (f>|), (|66|), and fl68|) run over 4 flavors of quarks 
(u, d, c, s), all assumed massless. At this energy, we do not include a b quark contribution in 
our calculation. For yfs = 58 GeV and 91 GeV, we use 5 flavors, again assuming all massless 
in the short- distance hard scattering cross sections. At these higher energies, non-zero mass 
effects for the c and b quarks are accommodated by our scale choice in the fragmentation 
functions, discussed below. 

The quark-to-photon fragmentation function that appears in Eq. (|63"D and (|68"D is ex- 
pressed as 



2.21 - 1.28.- + 1.29^ ^ , „ 

q l-1.63£n{l- z) v ; 



2vr 

xin (i4/vl) • (72) 



The gluon-to-photon fragmentation function in Eq. (B3) is 



z D g ^(z, /4) = °- 0243 (! " *) ^°- 9? • (73) 



These expressions for -D g ^ 7 and D g ^ y , taken from Ref. ||14j| , are used as a guideline for our 
estimates. The physical significance of scale /io is that the fragmentation function vanishes 
for energies less than /x . For g and for the u, d, s, and c quarks, we set [1q = Aq^ d , as in 



Ref. |Tj|. For the b quark we again use Eq. (|72|), but we replace /io by the mass of the quark, 
nib = 5 GeV; Db^(z, fj, F ) = for //p < m&. We set the fragmentation scale /ip equal to the 
renormalization scale \i for our inclusive cross sections. In the results presented below, we 
vary fj, to examine the sensitivity of the cross section to its choice. 

In presenting results, we divide our inclusive cross sections by an energy dependent cross 
section <7o that specifies the leading order total hadronic event rate at each value of y/s: 

Atts 
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(74) 



By doing so, we can observe what fraction of the total hadronic rate is represented by 
inclusive prompt photon production. 

In several figures to follow, we show the predicted behavior of the inclusive yield as a 
function of E 1 and 9 7 , as well as the breakdown of the total yield into contributions from 
various components. 

In Fig. |ll|, we present the inclusive yield as a function of E 1 at y/s = 91 GeV for two 

cLS cL 



values of the scattering angle 9 y , 45° and 90°. The same results are displayed in Fig. [12 
function of scattering angle 9 1 for two choices of E T In both Figs. |Tl| and 0, we set renor- 
malization/fragmentation scale fi = E T Dependence of the cross sections on \x is examined 



in Fig. [L3| at fixed E 1 . The patterns evident in Figs. pTTHT3| at y/s = 91 GeV are repeated with 



subtle differences in Figs. [T4"HT6| at y/s = 58 GeV, appropriate for experiments at TRISTAN, 
and in Figs. at y/s = 10 GeV, applicable for studies at CESR/CLEO. In Fig. [20] we 

compare predictions at the three energies by showing the cross section a^ l da / dx^dVL^ as a 
function of the scaling variable x 1 = 2E 1 /y/s. 



Evident in Figs. [TT] pT9] is the dominance of the lowest-order contribution to the inclusive 
yield, Eq. fl63|), at all values of y/s, except at small values of E^/y/s or at small values 
of n where the O (a em ) "direct" contribution, Eq. fl53|), becomes larger. Following the 
lowest-order contribution in importance at modest values of E 1 / y/s or of \x is the O (a em ) 
direct contribution. The direct contribution falls away more rapidly with increasing E 1 or 
\i than the O (a s ) quark-to-photon fragmentation term, Eq. (|68|). Therefore, at large values 
of Ey/y/s or fi, it is the O (a s ) fragmentation term, that is secondary in importance to 
the lowest-order term. The gluon-to-photon fragmentation contribution, Eq. fl6"5|) plays an 
insignificant role except at very small E T 

In Figs. [12], [15|, and [18], we examine the predicted # 7 dependence of our cross sections. 
These figures, presented with a linear scale, show perhaps more clearly the importance of 
the roles of the O (a em ) direct and O (a s ) fragmentation contributions. The lowest-order 
contribution, Eq. (|B"BP, is proportional to (1 + cos 2 # 7 ). However, there are significant sin 2 ^ 7 
components in the next-to-leading order direct term, Eq. (j64|) , and the next-to-leading order 
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fragmentation terms, Eqs. ( j66|) and (|68|). The net result is that the predicted total yield in 
Figs. [12], |T3j, and [18] is not proportional to (1 + cos 2 # 7 ). As illustrated in the figures, the 
deviation of the total yield from the (1 + cos 2 # 7 ) form becomes greater at smaller values of 
E 1 . (The results shown in Figs. [L^, [L5|, and [18] all pertain to the scale choice \i = E 7 .) One 
lesson from this examination of dependence of # 7 is that it is inappropriate and potentially 
misleading to assume that the functional form (1 + cos 2 # 7 ) describes the data when attempts 
are made to correct distributions in the region of small # 7 (where initial state bremsstrahlung 
overwhelms the final state radiation in which one is interested). 



Dependence on the renormalization/factorization scale \i in Figs. [13], |1| and [19| shows 
several interesting features. As is expected from the functional form of D q ^ (z, fi 2 ) in 
Eq. ([72]), the lowest-order contribution, Eq. (^), increases logarithmically as /x is increased. 
On the other hand, the in (s//x 2 ) dependent term in Eq. fl64|) causes a decrease of the O (a em ) 
direct contribution as \i is increased. Indeed, the (1 + cos 2 6* 7 ) part of the direct contribution 
becomes negative when sx 2 (1 — x 7 ) / jj? < 1. The physical cross section, represented as a 
solid line in Figs. |13|, [H| and [19], is of course always positive. 



An especially noteworthy feature of Figs. [L3[ [16], and |19| is that the total inclusive yield is 
nearly independent of /i, in spite of the strong variation with of its components. This in- 
dependence reflects the role of the fragmentation scale \i. It is introduced to separate "soft" 
and "hard" contributions into "fragmentation" and "direct" pieces. As the scale fi is in- 
creased, more of the cross section is necessarily factored into the fragmentation contribution, 
and vice versa, such that the sum remains nearly constant. 



In Fig. |20, we show the overall ^J~s dependence of our predictions. To facilitate com- 
parison, we present these results in terms of the "scaling" distribution a^da / dx^dVt^. The 
case of y/s = 10 GeV is somewhat special since we do not include a contribution from b 
quark fragmentation at this energy. Otherwise, the contribution of the lowest-order process, 
Eq. fl63|), decreases at fixed x 1 as a/s is increased. This decrease is explained easily. In com- 
puting da/dxrydfl-f, we multiply Eqs. ( |B"9"P and (|T2"P, obtaining a charge weighting factor of 
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egF^ c \s), whereas in computing the denominator ctq, the factor is F^ c (s). Owing to the val- 
ues of the v and a couplings in Table 1, the up- type quark contribution to F^ c (s) decreases 
as y/s increases, and the down- type contribution increases. The 0(a em ) direct contribution 
to gq 1 da I 'dXydQj decreases at fixed X*y clS yS is increased from 10 to 91 GeV. Again, the 
explanation may be found in the energy dependence of the ratio J2 q e2 q Fq C { s )l J2 q F^ c (s). 
Taken together these statements explain the energy dependence displayed in Fig. |(|. 

As remarked earlier, the particular expressions we chose for the fragmentation functions 
are not meant to be anything but illustrative expressions. It would be very valuable if 
these non-perturbative functions could be determined directly from data. Dominance of the 
q — > 7 fragmentation contribution in Figs. |TTHl9l demonstrates the important role data from 



e + e~ — > 7X may play in the extraction of D q ^ (z, fi 2 ) and study of its properties. However, 
as mentioned in the Introduction, an important limitation of high energy investigations is 
that photons are observed and cross sections are measured reliably only when the photons 
are relatively isolated. Since fragmentation is a process in which photons are part of quark, 
antiquark, and gluon jets, isolation reduces the contribution from fragmentation terms. In a 
forthcoming paper || , we will examine in detail the behavior of the isolated prompt photon 
cross section. 

In this paper we have presented a unified treatment of inclusive prompt photon pro- 
duction in hadronic final states in e + e~ annihilation. We have computed analytically the 
direct photon contribution through O (a em ) and the quark-to-photon and gluon-to-photon 
fragmentation terms through O (a s ). We presented the full angular dependence of the cross 
section, separated into transverse (1 + cos 2 ^) and longitudinal components. 
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APPENDIX A: TWO AND THREE PARTICLE PHASE SPACE 



In this Appendix, we express two- and three-particle final state phase space dPS^ and 
dPS^ in n dimensions in terms of the variables necessary for our calculation. We work out 
first the specific case of e + e~ — > qq. The four-vector momenta of q and q are p q and p q . 

The two particle phase space element in n = 4 — 2e dimensions is 

d n ~ l p q d n ^ P g 



dPS® 



.(2n) n 5^(q-p q -p q ). 



(Al) 



(2Tr) n - 1 2E q {2n) n ~ 1 2E q 
In the center of mass frame, p q = —p q and E q = E q . Eliminating the d n ~ l p q integration, we 



obtain 



dPS® 



1 

2^ 



n-2 



E n q ~ A dE q dd sin"" 3 9dtt n - 3 (p q ) 5 (e 9 - ^) • (A2) 



Since the square of the invariant matrix element, Eq. ( |17D of the text, depends on 6 but not 
on other angles, we may perform the integration over dQ n _ 3 ; 



n 



n— 3 



2tt 7T"7r(l - e). 



(A3) 



We derive 



dPS &> = i 1 d3p « 



2 (2tt) 3 E q 



An 



{s/A)sm 2 6j r(l-e) 



with x q = 2E q /\j .s. 



2tt 5(x q - 1) 
s x n 



(A4) 



For the three particle final state e + e~ —>■ qq^y, we label the four-vector momenta of q, q, 
and 7 as Pi,p%, and £. The invariant matrix element of interest to us, as defined in Eqs. ( p0[) 
and (p4|), depends explicitly on the inner products p\ ■ £,p% ■ £, and p\ - pi as well as on p\ ■ k, 
P2 ■ k, and I ■ k where k, defined in Eq. (|7]), is the difference k = k e + — k e - of the four- 
momenta of the initial e + and e~ . However, all these inner products are not independent. 
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Using momentum conservation and the fact that the momenta I and k are observables, one 
may show that p\ ■ I and p 2 ■ k are the only independent invariants. Note that it is completely 
equivalent to choose P2 instead of p\. 

For general orientation, it is useful to begin in n = 4 dimensions to establish the 
angular variables of integration we would use in that case, before generalizing to n di- 
mensions. In the overall e + e~ center of mass frame, we imagine a coordinate system with 
the 7 defining the z axis, vector k lying in the (x, z) plane, and vector f>\ generally having 
non-zero x, y and z components. 

k = \k\ (sia9 y , O,cos0 7 ) ; (A5) 

Pi = \Pi I ( s i n 0i7 cos 0; s i n sm 0; cos $17) j (A6) 
Pi ■ k = —pi ■ k 

— ~ \Pi\ \k\ (sin # 7 sin 6> l7 cos + cos^ cos 6 l l7 ) . (A7) 

The four-dimensional example shows that only the components of p\ in the £, k plane 
contribute to p ■ k. We use 9 X to denote the n-dimensional generalization of the four- 
dimensional azimuthal angular variable (ft, and we will express dPS^ in n-dimensions in 
terms of integrations over #i 7 and 9 X . 

Three particle phase space in n-dimensions is 

dpsm - Wr^k Wr=k {2 * n " Pl " P2 " iy (A8) 

Using the 6^ n ' function to eliminate the integrations over p 2 , we obtain 

dPS® = Pl , , -2tt5(^-E^-E 1 -E 2 ). (A9) 

Since we are interested ultimately in the invariant cross section E 1 da/d 3 £, we rewrite 
Eq. Q5g) as follows. 

cT"^ 1 



K™~ 2 d£ 7 dfl 7 sin n - 3 ^ 7 dfi n _ 3 (£) 



7 

"7 



(27r)"- 1 2 J B 7 (2tt)"- 1 2 j E; 7 

-7 — ^ — 7 — s — 7 (£ 2 sm 2 #J 2 — _^v_^_ A1Q 



2 (2tt) 3 £ 7 (2tt) 
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We take angle 9 1 to be the polar angle of the 7 with respect to the e + e~ collision axis in 
the overall center of mass frame. Since the square of the matrix element does not depend 
on 7 , we can integrate over dfl n -3(£) and c?0 7 independently. Using f2 n _ 3 from Eq. (|A3|), 
and / d(j) 1 = 2n, we reexpress Eq. ( |A10| ) as 



d n ~H 1 1 d 3 £ 
(2tt)™- 1 2E 7 ~ 2 (2tt) 3 eZ, 



47T 



EHirfeJ r(i-e) 



(All) 



We write cf™ Vi in Eq. ( |A9|) as 



with 



tT"^! = E?- 2 dE 1 dQ n - 2 ( Pl ) 

= El' 2 dE x dQ Xl sin"- 3 ^ l7 dQ n - 3 ( Pl ), 

dO n _ 3 (pi) = sin"" 4 ^ rffi n _ 4 (pi) 

= dcos6 x (1 — cos 2 ^)^ - <ifi n _ 4 (pi) 



(A12) 



(A13) 



Since only the components of p\ in the £, k plane contribute to pi ■ k, as shown in Eq. ( |A7| ), 
all angular variables on which the invariant matrix element depends are displayed explicitly 
in Eqs. ( |A12[ ) and ( |A13| ). We may therefore integrate dQ n -4(pi) in Eq. ( |A13| ) to obtain 



£4-4 (Pi 



2"- 4 7r^r(^) 



T(n - 4) 

In this frame, E2 in Eq. (|A9|) can be expressed as 

E 2 2 = (P2) 2 = (pi + t? = El + E 2 + 2E 1 E 1 cos9 



(A14) 



it 



(A15) 



Using Eq. ( |A15|) , we can replace the integration over <icos6 l i 7 in Eq. (|A12|) by an integration 
over dE 2 ; for fixed E\, 



EiE^dcosOi^ = E2 dE2 



(A16) 



Substituting into Eq. 



we derive 
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dPS (3) 



1 1 d 3 £ 



47T 



2 (2tt) 3 E 1 [{E^ sin 2 6 J r (!- e ) 



1 1 



4vr 2 



X 



4 (2tt) 2 V^isin 2 ^ 



^ fi n _ 4 ( Pl )(l-cos 2 ^)- £ - 1/2 rfcos^ 



x — 5{Vs - E y — Ei — E 2 ) dE l dE 2 



E, 



(A17) 



It is easy to verify that Eq. (|A17 ) reduces to the familiar form in four-dimensions when 
e -> 0. 

We introduce new dimensionless variables related to the singularity structure of the 
invariant matrix elements, given in Eqs. fl3"Dp and (3"1|): 

2pi • £>2 



1/12 

2/2^ 
X\ 

x 2 



2pi-£ 
q 2 

q 2 
2gi ■ g 

g 2 

2p 2 • g 

q 2 

21 -q 



(Al8a) 
(A18b) 
(A18c) 
(A18d) 
(Al8e) 
(A18f) 



We observe that yie = 1 — ^2; 2/2^ = 1 — and ?/i2 = 1 — x T In the center of mass frame, 
q = (\/s, 0), we have 

2E 2 



2K 



Xi 



and x 2 



(A19) 



s s 
dEi dE 2 = — dx\ dx 2 = - dyu dy 2 £ 



(A20) 



5(y/s - E q - E q - E y ) = —i=d(l - y u - yn - yu) 

v s 



After some algebra, one may verify that 



(A21) 



El E 2 sin 2 9 



It 



yu yu 2/i2- 



(A22) 
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Substituting Eqs. (g2g )-( g2g) into Eq. flATT] ), we derive 



dPS {3) 



1 1 dH 



2 (2tt) 3 



47T 



(s/4)sin 2 # 7 / r(l-e) 



x 



\ j^y 2 (^f) ^u-a(Pi) (l-cos^^-^dcos^ 
%u cfc/2<? 8{l - yu - y 2 t - 2/12) 



2/2£ 2/12)' 



For reference, we record that 



dcos9 r (l — cos 2 6L) e 



r(|)r(|-e) 
r(i-c) 



(A23) 



(A24) 



dcos8 r cos 6L (1 — cos 9„) 



r(l)r(|- e ) 



r(2-e) ' 

Dividing Eq. ( |A25|) by Eq. ( |A24|) , we define the average of cos 2 8 X in n-dimensions as 

1 1 



(A25) 



(cos 9 X ) 



2 1-e' 



(A26) 



When e — > 0, 9 X becomes the azimuthal angular 0, and Eq. ( |A26| ) is consistent with the 
4-dimensional result 



1 r 27T 



cos 0) 



cos 



(A27) 



27rJo 2 

Equation (|A23|) is written in a form for which the photon with momentum I is the 
observed particle, with integrations done over the momenta of other final state partons. 
When considering e + e~ — > qqg with q (or q) fragmenting into the observed 7, we require 
instead E q da/d 3 p q . It is useful, therefore, to reexpress Eq. (|A23|) in a form that manifests 
the symmetry of phase space among all three final state particles. 

Using yi^ = 1 — x 1 , we introduce the identity 



1 = dy 12 5 (x 7 - (1 - 2/12)) 



(A28) 



Inserting this identity into Eq. ( [A23j ), we obtain a more symmetric form 



33 



4 



2tt 



rn 



dyu dy 2 e dy 12 , , 
x 0(1 - y u - Vii - 2/12) ■ 



(A29) 



Vu Vu Ui2 

We can use Eq. (|A29D to derive an expression for three-particle phase-space suitable for 
calculating E q da/d 3 p q when the quark from e + e~ — > qqg is the fragmenting parton. We 
label (q,q,g) as (1,2,3). We let the gluon replace the photon 7 in Eq. (|A29|) , and we then 
switch the labels of the g and the observed quark. We derive 



dPS^ = I ' 



d 3 px 



x 



2 (2tt) 3 E l 
s 



An 



[si A) sin 2 61 



r 1 



2n S(xi 
s 



2/23)) 



Xi 



1 \ 2 /An 



dQ n _ 3 (p 3 ) 



2ttJ \ s J T(l - e) a n _ 3 ( P3 ) 



dy\2 dy 13 dy 23 
x — : — 5(1 - y 12 - 2/13 - J/23) • 



(A30) 



2/12 2/13 2/23 

The first line of Eq. ( \K30j ) is identical to dPS {2 \ Eq. (0). By switching 1 and 2 in Eq. ([MUD , 
one can get a form for the phase space suitable for calculating the hard-scattering cross 
section when the antiquark is the fragmenting parton. 
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FIGURES 

FIG. 1. Illustration of e + e~ — > in an m parton state: e + e~ — ► followed by fragmentation 

FIG. 2. Contribution to the hard-scattering cross section for e + e~ — > cX with the fragmenta- 
tion of c — > 7X. Here c denotes an intermediate photon or a gluon or a quark of any flavor. 



FIG. 3. Lowest order, O (a° m a°), photon production through quark fragmentation. 

FIG. 4. Feynman diagrams for e + e~ — > ~/qq. 

FIG. 5. Order a s Feynman diagrams for e + e~ — > that contribute to e + e~ — > 7X via <? ^ 7 
fragmentation. 

FIG. 6. Contributions to the 0{a s ) cross section cr^+ e -^ qX ; (a) real gluon emission diagrams 
(e + e~ — > (/gg), (b) virtual gluon exchange diagrams that interfere with the lowest order tree 
diagram. 

FIG. 7. Order 0(al rn ) contribution to the hadronic tensor H^ v . 

FIG. 8. Center of mass coordinate axes of an e + e _ collision with the z-axis being the direction 
of the observed photon. 



FIG. 9. Feynman diagram for D^(z). 



FIG. 10. Feynman diagrams for D\tl q {z). 
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FIG. 11. Normalized invariant cross section for the inclusive preocess e + e~ — > 7X at 
y/s = 91 GeV shown as a function of the photon energy E~ for two values of the photon scat- 
tering angles: a) # 7 = 45°, and b) # 7 = 90°. Displayed are the total result and the four separate 
contributions from lowest-order fragmentation ( "oth-frag" ) , 0(a em ) direct production, and the 
0(a s ) quark and gluon fragmentation contributions. The renormalization/fragmentation scale 
fi = Ej. 

FIG. 12. Normalized invariant cross section for the inclusive preocess e + e~ — > 7X at 
y/s = 91 GeV shown as a function of the photon scattering angle # 7 for two values of the photon 
energy: a) Ey = 15 GeV, and b) Ey = 30 GeV. The total result and the four separate component 
pieces are displayed. The renormalization/fragmentation scale = E~. 

FIG. 13. Renormalization/factorization scale dependence of the normalized invariant cross sec- 
tion for the inclusive process e + e~ — * 7V at y/s = 91 GeV for two values of the photon scattering 
angle: a) # 7 = 45°, and b) By = 90°. The photon energy Ey = 15 GeV. The total result shows little 
\i dependence, whereas the component contributions display considerable compensating variation 
with jj,. 



FIG. 14. Photon energy dependence, as in Fig. 11, but for center-of-mass energy y/s = 58 GeV 



FIG. 15. Photon scattering angle dependence, as in Fig. [T^, but for center-of-mass energy 
y/s = 58 GeV, and photon energies a) Ey = 10 GeV, and b) E y = 20 GeV. 



FIG. 16. Renormalization/factorization scale dependence, as in Fig. 13, but for center-of-mass 
energy y/s = 58 GeV, and photon energy E^ = 10 GeV. 



FIG. 17. Photon energy dependence, as in Figs. |ll| and 14, but for center-of-mass energy 
y/s = 10 GeV. 

FIG. 18. Photon scattering angle dependence, as in Figs. ^ and but for center-of-mass 
energy y/s = 10 GeV, and photon energies a) = 1.5 GeV, and b) Ey = 3 GeV. 
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FIG. 19. Renormalization/factorization scale dependence, as in Figs. 13 and 16, but for cen- 
ter-of-mass energy y/s = 10 GeV, and photon energy Ej = 2 GeV. 

FIG. 20. Normalized invariant cross section for the inclusive preocess e + e~ — > 7X expressed 
in terms of the scaling variable x 7 = 2E^/y/s. Results are presented for two values of the photon 
scattering angle: a) 9 1 = 45°, and b) # 7 = 90°. Shown are curves for three center-of-mass energies: 
y/s =10, 58, and 91 GeV. 
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TABLES 



TABLE I. Electroweak V-A coupling constants 




(-1 + 4 sin 2 W ) /(2sin 


26 w ) 


a e 


1/(2 sin 26 w ) 




v q 


(7| - 2e q sin 2 6 W ) / (sin 


26 '„,) 


a q 


-III sin 26 w 




TABLE II. Isospin and fractional charges for quarks 




4 


e q 


u, c, t 


1/2 


2/3 


d, s, b 


-1/2 


-1/3 
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Fig.lB-(b) 
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Fig.l9-(a) 
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Fig.l9-(b) 
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